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We consider the effect of future dark energy experiments on "Albrecht-Skordis" (AS) models of 
scalar field dark energy using the Monte-Carlo Markov chain method. We deal with the issues 
of parameterization of these models, and have included spatial curvature as a parameter, finding 
it to be important. We use the Dark Energy Task Force (DETF) simulated data to represent 
future experiments and report our results in the form of likelihood contours in the chosen parameter 
space. Simulated data is produced for cases where the background cosmology has a cosmological 
constant, as well as cases where the dark energy is provided by the AS model. The latter helps us 
demonstrate the power of DETF Stage 4 data in the context of this specific model. Though the AS 
model can produce equations of state functions very different from what is possible with the wo — w a 
parametrization used by the DETF, our results are consistent with those reported by the DETF. 

PACS numbers: 



I. INTRODUCTION 



Current astronomical observations are best fit by cos- 
mologies containing a significant density of energy with 
negative pressure. This has been dubbed dark energy, 
and its exact nature remains a mystery. Though the 
properties of the dark energy are poorly constrained by 
current data, a number of proposed experiments will 
have the reach to probe them. Most assessments of pro- 
posed experiments use abstract parameters such as the 
11 wq — Wa' 1 parametrization of the dark energy equation of 
state used by the Dark Energy Task Force (DETF) 1], 
or generalized forms of scalar field potentials 0,0, 4|. 
In order to fully evaluate future experiments it is use- 
ful to understand the impact and discriminating power 
they will have on specific proposed models of dark en- 
ergy. In this paper, we will focus on one scalar field 
(or "quintessence") model of dark energy, the so-called 
"Albrecht-Skordis" (AS) model Hi]. 

The AS model is interesting for a number of reasons. 
Fields of the form needed for the AS model can arise as 
radius moduli of curled extra dimensions 0, and, while 
these moduli may be treated as scalar fields in gravi- 
tational calculations, they are not subject to the quan- 
tum corrections that can create serious problems for most 
quintessence models [1, i] . The phenomenology of these 
models is also interesting because, unlike the majority of 
dark energy models, the AS quintessence field can con- 
tribute a significant fraction of the total cosmic energy 
density throughout the history of the universe. 

Central to this work are the projected data sets (or 
"data models") created by the DETF 1]. We use the 
supernova, weak lensing, baryon oscillation, and Planck 
data sets (though not the cluster data sets, for techni- 
cal reasons similar to those outlined in 10]), to show 
how these data sets would impact the range of possi- 
ble parameters for quintessence models. To this end, we 
use Monte-Carlo-Markov chains (MCMC) on these data 
sets to analyze dark energy and cosmological parameters, 



rather than the fisher matrix methods used by the DETF. 
In addition to using "standard" simulated data based on 
a universe with a cosmological constant, we create pro- 
jected data sets based on an AS model universe, and run 
the MCMC around these to highlight the discriminating 
power of the projected data sets. This paper is one in a 
series of papers that consider different scalar field models 
using similar methods [ll|. (Technical information about 
our data models and MCMC methods will be presented 
in an appendix of another paper in the series [111]). 

We show that the stage-by-stage improvement in pa- 
rameter constraints for the Albrecht-Skordis model is 
similar to the relative improvement in the wq — w a con- 
straints seen by the DETF. We do the same for data sets 
generated around a specific AS model which shows small 
deviations from w(a) = —1 in the present epoch. This 
AS model is chosen to illustrate a case where if the real 
universe is described by an AS model, pure cosmological 
constant dark energy can be ruled out by a large mar- 
gin using good Stage 4 data. We also note some features 
in the AS parameter contours when fitting to cosmologi- 
cal constant based data that reflect the fact that the AS 
models we consider can only ever duplicate a cosmologi- 
cal constant in an approximate manner. 



II. PARAMETER SPACE OF THE 
ALBRECHT-SKORDIS MODEL 

Scalar field models of dark energy, or quintessence, are 
considered in the framework of an FRW cosmology, with 
and equation of motion of 



2-(j) + a -— = 
a od> 
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where a is the scale factor, <f> is the homogeneous scalar 
field, and V its potential. The background is assumed to 
be the homogeneous Friedmann equation 
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The equation of state parameter of the scalar field is 



to 



P<t> 
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where p<j, and are the pressure and density of the scalar 
field. Non-relativistic matter has w — 0, radiation has 
w = 1/3, and a cosmological constant has w = — 1; a 
homogeneous scalar field can have any behavior in this 
range, depending on the actual time evolution <p(t). The 
Albrecht-Skordis model[|[ postulates a scalar field with 
an "exponential-with-pre-factor" potential, 
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One of the attractions of this model is that realistic cos- 
mologies can be achieved for cases where the parameters 
in this potential are all roughly of order unity in Planck 
units. Specifically, with such parameter choices the po- 
tential can have a local minimum with a height consis- 
tent with the dark energy density observed today. It is 
assumed that the field starts high up on the potential. 
The field will then rapidly approach an attractor solu- 
tion where the field mimics the equation of state of the 
dominant energy. It continues this tracking behavior un- 
til the field reaches the (approximately) quadratic shaped 
local minimum, at which point it begins a damped oscil- 
lation around that minimum, giving its equation of state 
a wave like variation as the field becomes the dominant 
energy. Figure [1] shows w(a) throughout the history of 
the universe for a typical AS model. 

The potential of the AS model has in principle four de- 
grees of freedom important to cosmology, most of which 
are complicated combinations of the various parameters. 
Figure [2] illustrates some of the parameter dependences 
of the potential. Perhaps the most important degree of 
freedom, because it determines the scalar field fraction 
of the total energy density during the tracking behavior, 
is the logarithmic slope, (InV) 1 or this is almost en- 
tirely controlled by A. The height of the local minimum 
sets the dark energy density today. The curvature of the 
potiential at its minimum is related to the height of the 
local maximum. The height of the local maximum deter- 
mines whether the field stops in the minimum or rolls off 
to infinity, and the curvature at the minimum sets the 
frequency of the late time oscillations. We have found 
that the oscillation frequency is not measurable by any of 
the simulated data sets we considered. Also, the param- 
eter space for models that roll off to infinity but still give 
realistic cosmic acceleration is a only an exponentially 
thin and extremely finely tuned region We choose 
to ignore this exotic behavior and focus completely on 
parameters for which the field (classically) stays in the 
local minimum. 

It will be helpful to re-express Eqn. 0] in the following 
parameters: 
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FIG. 1: The equation of state w for the scalar field of a typical 
AS model. Note that here the a scale is logarithmic in order 
to show behavior on all time scales. When the background 
energy density drops to a level approaching the initial energy 
of the field during radiation domination, the energy density 
goes through a transient before mimicking the equation of 
state of radiation (w = 1/3). At matter domination, the field 
again goes through a transient, but the field approaches the 
local minimum in the potential and becomes the dominant 
form of energy (with w — ► — 1) before it can stabilize with a 
mater-like equation of state (w = 0). 
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FIG. 2: The effects of the original parameters (from Eqn. [4]) 
on the potential minimum. All measures are in Planck units. 
Note that very small changes in A or B cause substantial 
changes in the height of the minimum. 



Because the number of parameters is greater than the 
number of parameters measured by our simulated data 
the AS model has a number of degenerate directions. 
Along these degenerate directions, the parameters can 
change without significantly affecting the observables. 
The most apparent of these directions are the ways that 
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FIG. 3: Under the new parameterization, Vo maintains it's 
previous effects, while the new parameter f3 becomes the pri- 
mary parameter for the AS model. (3 sets the exponential 
slope as well as the minimum curvature. 
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FIG. 4: This figure illustrates how the primary parameter /3 
affects the w(a) behavior of the AS model. The AS model- 
based data was generated using (3 = 80. 



the height of the local minimum of the potential can be 
changed. This height is exponentially sensitive to the 
product of the exponential factor A, and /?, which pri- 
marily controls the location of the minimum. This de- 
generacy can be stabilized by fixing the product of A and 
/3, so that A = There is also a dependence of the 

,2 

minimum on 5 and — . To work around these difficulties, 
x _ ' 

we fix 5, and vary x as a function of (3 so that the overall 
constant Vo is the only parameter that adjusts the height 
of the minimum. This necessarily excludes the portion of 
the parameter space where there is no minimum (with es- 
sentially no loss of generality, as discussed above). Also, 
these parameter constraints do not allow the oscillation 
frequency near the local minimum to vary independently 
of the other parameters. This feature is required to suc- 
cessfully run the MCMC calculations since none of the 
data sets are sensitive to the oscillation frequency. Figure 
[3] shows how j3 affects the new parameterization. 

Aside from degeneracy in the specifics of the potential, 
there are broad patches of parameter space that are, for 
most data sets, not detectably different from a cosmolog- 
ical constant. This happens when the exponential slope 
is steep, and the field slides very early down to the local 
minimum and oscillates in a shallow and rapid manner, 
resulting in an equation of state very nearly —1 at the 
times of interest. 

The starting value of the scalar field is mostly unim- 
portant in the AS model, but it should be noted that 
there are values of (3 and the initial (f> that can cause 
problems. For large values of j3, if the scalar field starts 
at a point higher on its potential than the radiation den- 
sity during nucleosynthesis, then the scalar field will form 
a significant fraction of the energy density at that epoch. 
Such regions of parameter space would be ruled out by 



observations of primordial element abundances, though 
our algorithm does not contain such considerations. 

Other related problems with the algorithm specific to 
the AS model include the last scattering surface. The 
algorithm does not calculate the size of perturbations at 
the time of last scattering, but takes them as an input, 
so the effect of the scalar field on that size is not taken 
into account by the algorithm even when the scalar field 
energy density is significant. With respect to the growth 
of linear density perturbations, the algorithm does not 
use the scalar field in its calculation until a redshift of ten. 
The effect of these issues is that some constraining power 
has been ignored. Albrecht and Skordis investigated this 
issue @ to some degree, but while we do not believe this 
affects the conclusions of this paper, this matter remains 
interesting and worthy of further exploration. 

In our parameterization, the amplitude of the oscilla- 
tions in the equation of state at late times is primarily 
determined by the value of (3. For small values of (3 the 
oscillation amplitudes are very small, but become large 
for values of (3 approaching 100. This somewhat unintu- 
itive behavior can be explained by noting that, when the 
dominant form of energy is scaling as a~ n , p$j p c = jj 
(where p c is the critical density) [13J] . Thus, the smaller 
(3 is (and the larger A is), the further the field energy is 
below the matter density during mater domination, and 
the more Hubble friction will slow the field. The closer 
the scalar field energy density is to the matter density 
during matter domination, the more kinetic energy re- 
mains in the field when it becomes dominant, and thus 
free to oscillate. 
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FIG. 5: Likelihood contours for the AS model in f3 and Vo 
for Stage 2 (top) and Stage 3 photometric, optimistic (bot- 
tom), for data generated from a cosmological constant model. 
Though constants in Vo improve greatly, the AS model has 
little observable difference from a cosmological constant for 
values of f3 below 50. 



III. CONSTRAINING THE AS MODEL 
AROUND COSMOLOGICAL CONSTANT 
MODEL DATA SETS 



Having chosen to explore the scalar field parameter 
space in terms of the parameters j3, Vo, and the initial 
value of the scalar field (with the other parameters con- 
strained) , we then performed an MCMC analysis on sim- 
ulated data sets generated around a cosmology with a 
cosmological constant. The technical details of this work 
are presented in the appendix of [ll[ . These data sets are 
consistent with those used by the DETF, using simulated 
Stage 2, Stage 3, and Stage 4 supernova, weak lensing, 
baryon oscillation, and CMB observations. We did not 
use the cluster data because of the difficulty of adapt- 
ing the DETF construction to a quintessence cosmology, 
nor have we included possible improvements from cross- 



FIG. 6: Likelihood contours for the AS model in /3 and Vo 
for Stage 4 space (top), and Stage 4 ground LST (bottom), 
both optimistic cases, for data generated from a cosmological 
constant model. Though constants in Vo improve greatly, the 
AS model has little observable difference from a cosmologi- 
cal constant for values of /3 below around 50. The different 
regions represent 68.27%, 95.44% and 99.73% confidence re- 
gions. This convention is used for all contour plots in this 
paper. 



correlations [14j, [l5|. In the DETF language Stage 2 
encompasses projects that are fully funded and under- 
way, Stage 3 models medium term, medium cost propos- 
als, and Stage 4 are the larger projects, such as a large 
ground-based survey or a new space telescope. 

Though there are three parameters controlling the 
state of the scalar field, the initial value of the scalar field 
is generally washed out by the tracking behavior, and the 
variance in height of the minimum is mostly a reflection 
of the uncertainty in the Hubble parameter. This leaves 
us with one important parameter, [3, that mostly controls 
the unique properties of the AS model. We include a Vb 
axis in the contour plots for convenience. One should 
take care, though, to remember that Vo is not a param- 
eter that can be said to have a significant effect on the 
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equation of state of the dark energy. 

Figure [5] shows the likelihood contours for DETF Stage 
2 and Stage 3 (optimistic) data. From the point of view 
of these data sets the behavior of the AS model can be 
thought of as trending toward a cosmological constant 
as (3 goes toward zero. While these plots show that the 
constraints on Vo, which corresponds to the dark energy 
density now, improve greatly with stage number, the crit- 
ical consideration here is (3. The constraints on f3 do not 
appear to improve in as dramatic a manner. This is be- 
cause there is a substantial range of (3 where the these 
data sets cannot distinguish the AS model from a cosmo- 
logical constant. As seen in Fig [4] at (3 — 50 already one 
must look back to a = .2 (z = 4) to see a large deviation 
from w = — 1. To some extent this feature will affect 
any attempt to evaluate the AS model using data from 
a universe with a cosmological constant. In that low (3 
region, the small oscillations of the field persist at late 
times, and the Stage 4 data sets may be strong enough 
to pick up on some effects of these oscillations. Figure [5] 
shows likelihood contours for DETF Stage 4 space and 
ground (optimistic) data. There is a kink in the Stage 
4 space contour plot which corresponds to a range of (3 
where the residual oscillations at late times peak in am- 
plitude. As the oscillations represent an energy density 
above the amount corresponding to the minimum of the 
potential, Vq must be lower to accommodate regions with 
larger late time amplitudes. 

Figure [7] shows contours in ujk — f3 space (tOk , the "cur- 
vature density" is a measure of the curvature of the uni- 
verse). There are well known degeneracies associated 
with determining curvature and the dark energy equa- 
tion of state simultaneously [16| . Here these degeneracies 
show up by allowing "best fits" to data from a flat uni- 
verse using non-zero curvature for certain values of (3, as 
can be seen in contours for the Stage 2 data. The Stage 
4 data set shown in the lower panel clamps down on this 
behavior considerably. 



IV. MCMC OF THE AS MODEL AROUND AS 
MODEL-BASED DATA SETS 

As a way of further exploring the power of the ad- 
vanced data sets, a sample AS model was chosen as a 
new fiducial model around which the data sets were gen- 
erated. This model was chosen to be marginally consis- 
tent with the cosmological constant-based Stage 2 data, 
but different enough from a cosmological constant that a 
cosmological constant would be strongly ruled out with 
Stage 4 data (the point being to illustrate the power of 
Stage 4). The fiducial model chosen here has a (3 = 80, 
as depicted in Fig. 2] Because the Stage 2 data allowed 
larger values of (3 when curvature was present, the sample 
AS model was chosen with non-zero spatial curvature to 
make it more consistent with a flat cosmological constant 
model for this value of (3. The results of the calculation, 
as seen in Fig. [5] and show that, for a universe de- 
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FIG. 7: This figure illustrates the confusion of the AS model 
with curvature. The Stage 2 set (top) allow significant de- 
partures from flatness, even when using data generated by a 
cosmological constant model with no curvature. The Stage 4 
(ground, optimistic) data set (bottom) are effective at con- 
straining this degeneracy. 



scribed by this specific model, the Stage 4 experiments 
will rule out a cosmological constant by several sigma. 



V. DISCUSSION AND CONCLUSIONS 

We have analyzed the impact of the DETF simulated 
data sets in the context of the Albrecht-Skordis scalar 
field model of dark energy. We find that the effect of 
the DETF data sets on the parameter space of the AS 
model is very similar to the DETF findings. There is sub- 
stantial improvement in the constraining power of each 
successive stage of experiments. We presented likelihood 
contour plots in AS parameter space for a few key com- 
binations of DETF simulated data. We used these plots 
to demonstrate the broad agreement with the DETF re- 
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FIG. 8: Likelihood contours for the AS model in /3 and Vb for 
Stage 2 (top) and Stage 3 photometric, optimistic (bottom), 
for data generated from the j3 — 80 AS model. 



FIG. 9: Likelihood contours for the AS model in /3 and Vb 
for Stage 4 space (top), and Stage 4 ground LST, both op- 
timistic (bottom) cases, for data generated from the selected 
AS model. 



suits and also to point out more subtle effects such as de- 
generacies between the dark energy equation of state and 
curvature and some peculiarities of the way the AS model 
can approximate a cosmological constant. In the course 
of this work we have produced similar contour plots for 
many more instances of the DETF simulated data (taking 
individual techniques separately for example, and using 
both the DETF optimistic and pessimistic data models). 
We find our overall conclusions and detailed points apply 
quite generally across the full range of DETF simulated 
data. We also find that as one moves to better experi- 
ments, the improved discriminating power will definitely 
rule out interesting portions of the AS model space and 
could completely rule out the cosmological constant if the 
AS model is correct. 

The DETF reported a figure of merit in terms of the 
inverse of the area inside of likelihood contours in their 
model space, wq and w a . This figure of merit showed an 
improvement of a factor of three going from Stage 2 to 
Stage 3, and a factor of ten going from Stage 2 to Stage 



4, assuming good Stage 3 and Stage 4 projects. In the 
parameter space we used for the AS model we have only 
one parameter that strongly affects the equation of state. 
The constraints on this parameter, f3, show improvement 
by roughly the square roots of these factors, appropri- 
ate for the reduced dimensionality in the AS parame- 
ter space. We also see some subtle differences between 
the constraining power of ground and space Stage 4 data 
models in that the other variable plotted, Vo, is somewhat 
more strongly constrained by the Stage 4 ground data 
set than for the Stage 4 space. Other scalar field models 
[TTI . \vf\ display the opposite behavior. We are currently 
investigating this effect further, which may lead to inter- 
esting insights into the complementarity of ground and 
space-based Stage 4 experiments. 

It is interesting to consider the relationship between 
this work and recent work by one of us and G. Bernstein 
There the same DETF data sets were studied in 
the context of an abstract dark energy model where the 



7 



equation of state was modeled by many more than the 
two parameters used by the DETF. One conclusion of 
[lol ] was that high quality data can make good measure- 
ments of significantly more than two equation of state 
parameters. We have seen in this paper that the main 
equation of state parameter of the AS model is con- 
strained by DETF data sets to a similar degree as the 
DETF w a — w a parameters, even though they describe 
very different functions w(a). As discussed in [ijj], we 
believe that this ability to constrain a wide variety of 
functions w(a) is another manifestation of the rich con- 
straining power demonstrated in [icl |. 

A good way to think about this may be to consider the 
expansion of the AS family of w(a) function as well as the 
wo — w a family of functions in terms of the independently 
measured orthogonal functions Wi (a) from [Ic| . The fact 
that we have seen here (and also in [ll|, 03) that a vari- 
ety of different w(a) functions can be constrained as well 
as the DETF u>o — w a functions appears to reflect the 
fact that fundamentally many more functions w(a) are 
measured than are contained in any one of these families 
alone. 

An upshot of this and the companion work in (lH Il7|) 



is that modeling the impact of future experiments using 
the DETF parameters appears to give a good indication 
of the impact on quintessence dark energy models with a 
similar number of parameters. An advantage of the meth- 
ods used here is that one can see explicitly how future 
data can constrain real dark energy models in a signifi- 
cant way, and can even eliminate some models entirely. 
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